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Abstract 

By matching 1/m? divergences in finite- volume two-point correlation functions of the scalar or 
pseudoscalar densities with those obtained in chiral perturbation theory, we derive a relation 
between the Dirac operator zero-mode eigenfunctions at fixed non-trivial topology and the 
low-energy constants of QCD. We investigate the feasibility of using this relation to extract 
the pion decay constant, by computing the zero-mode correlation functions on the lattice 
in the quenched approximation and comparing them with the corresponding expressions in 
quenched chiral perturbation theory. 
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1. Introduction 



It is well appreciated that a combination of lattice methods and chiral perturbation theory 
(%PT) can be an efficient tool for studying the low-energy properties of QCD close to the chiral 
limit. While xPT is the perfect book-keeping device for the non-trivial relations implied by 
chiral symmetry, the lattice can be used to determine the low-energy couplings of this theory, 
which encode the dynamics of the fundamental Lagrangian. 

The study of QCD on the lattice obviously requires a finite volume, and this might ap- 
pear problematic close to the chiral limit, since spontaneous chiral symmetry breaking does 
not take place in a finite volume. This is not the case, however, because is able to 

predict analytically the large finite-size effects expected in this regime, in terms of the same 
low-energy constants as appear in an infinite volume, in such a way that infinite-volume 
quantities can be obtained unambiguously from the finite- volume ones. The study of xPT 
in a finite volume and close to the chiral limit (in the so-called e-regime) was pioneered by 
Gasser and Leutwyler PQ-[1] a long time ago, but it is only recently that practical "mea- 
surements" of physical observables became feasible in lattice QCD 5 - 8 . This is thanks to 
the new formulations of lattice fermions, which preserve an exact chiral symmetry 
In this paper we will employ one of these formulations and invoke the specialized numerical 
techniques developed in ref. ^7j> which are needed for high-precision studies in the e-regime. 

It was found in ^Hj that in the e-regime, gauge field topology may play a very important 
role. In a given chiral regularization of QCD, averages can be defined in sectors of fixed 
topological index v |19j , and our assumption will be that standard ultraviolet renormalization 
also makes sense in such sectors. Although this is a non-trivial assumption in QCD, there is 
a well-defined prescription for how to compute analogous averages in xPT. 

It then turns out that close to the chiral limit, many observables depend quite strongly 
on the topology. In particular, for v ^ 0, two-point functions of the scalar and pseudoscalar 
densities have poles in the quark mass squared, with residues given by correlation functions of 
Dirac operator zero-mode eigenfunctions. In the e-regime of xPT the same poles appear, with 
residues that are calculable functions of the low-energy constants. Requiring the residues in 
the fundamental and effective theories to be the same yields non-trivial relations. 

To be more specific, at leading order in xPT the correlators mentioned are constants 
depending only on v and the volume, but at next-to-leading order (NLO) one obtains a space- 
time-dependent function, which also involves the pseudoscalar decay constant F. Therefore, 
F can be determined by monitoring the amplitude of the time dependence. A nice feature 
of this procedure is that it does not require knowledge of renormalization factors since we 
employ a regularization that preserves the chiral symmetry. 

Given that F 2 appears first at the NLO, 0(1/ F 2 ), and that the convergence of xPT at 
realistic (not very large) volumes is not a priori guaranteed to be rapid, it is one of the 
purposes of this paper to present the results of the calculation up to the next-to-next-to- 
leading order (NNLO), 0(1/F 4 ). According to our conventions as detailed in Appendix A 
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these correspond to the relative orders 0(e 4 ),0(e 8 ) in the e-expansion, respectively. 

We study, furthermore, the feasibility of using this relation to extract F 2 from the zero- 
mode wave functions computed on the lattice, in the quenched approximation. Thus, pre- 
dictions for the quenched version of xPT (QxPT) (whose theoretical status is unfortunately 
rather questionable, see Sec. 13.11 and, e.g., ref. [20]) are also presented, at the same order. 

The paper is organized as follows. In Sec. [21 we derive the relation alluded to above and 
present the results of the calculation of the pseudoscalar density correlator in full xPT. In 
Sec. El we obtain the same results in the quenched approximation and compare them with 
a numerical determination of the zero-mode eigenfunctions in lattice QCD, using overlap 
fermions. We conclude in Sec. |1J and collect various details of the NNLO computations in 
three Appendices. 

2. Pseudoscalar correlator in QCD and in %PT 
2.1. The fundamental theory 

In this paper we are concerned with QCD in a finite volume V = T x L 3 , with periodic 
boundary conditions in all directions. Our conventions for the Dirac matrices are such that 
7/1 = 7(U> {ifiilu} = 2<5 M „, 75 = 7o7i7273 ; so that the (unquenched) Euclidean continuum 
quark Lagrangian formally reads 

£ E = ^D^ + M)i> , (2.1) 

where M is the mass matrix. For simplicity, we take M to be diagonal and degenerate, 
M = diag(m, ...,m). The number of dynamical flavours appearing in ip is denoted by Nf. 

In the following we will restrict our attention to correlation functions of the scalar and 
pseudoscalar densities, 

S 1 = tpT 1 ^ , V 1 = ^75T 7 V , (2.2) 

involving N v valence quarks; in the unquenched theory, N v = Nf. The iV v x N v valence 
flavour basis is generated by 

T I = {T Q ,T a }, T° = I Nv , a = l,...,iV v 2 -l, (2.3) 

where In v is the iV v x iV v identity matrix, and the traceless T a are assumed to be normalized 
so that 

Tr [T a T b ] = -5 ab . (2.4) 

Our analysis is based on the assumption that correlation functions at fixed topology, e.g. 
the two-point correlators of pseudoscalar densities, 

C I u J (x-y) = (v I (x)V J (y)) u , (2.5) 
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have a well-defined meaning in the continuum limit at non-zero physical distances. Although 
plausible, this is a non-trivial dynamical issue and to pose precise questions we must introduce 
an ultraviolet regularization. 

We here adopt the lattice regularization with a massless Dirac operator D obeying the 
Ginsparg- Wilson (GW) relation, since it preserves an exact chiral symmetry. The topological 
index assigned to a configuration then is v = n + — n_, where n + (n_) are the numbers of 
zero-modes of D with positive (negative) chirality. Correlation functions such as Eq. 1)2,5)) 
are now well defined at fixed cutoff an d the question is whether, in any given sector 
of index z^, they have a continuum limit independent of the particular choice of D . Our 
working hypothesis is that this is indeed the case; some recent numerical evidence (in the 
quenched approximation) consistent with this scenario can be found, e.g. in refs. [51 121j. 

By employing the spectral representation of the quark propagator, it is clear that the 
correlator in Eq. 1)2. 5 J) contains a pole in m 2 , due to the exact zero modes. Its residue is 

lim (mV) 2 Ci J {x) = Tr [T 1 ^] A u (x) + Tr [T J ]Tr [T J ] A v {x) , (2.6) 

where 

Au(x-y) = (J2 v}(x)vi(x)vl(y)vj(y)) v , (2.7) 

and the sums are over the set of \u\ zero modes Vi of the Dirac operator, Dvi = V i G K,, 
which have definite chirality and are assumed to be normalized so that Jd 4 x vj(x)vi(x) = 
V. Eq. 1)2. 7|) corresponds to a "connected" contraction of the quark lines, Eq. 1)2.8)) to a 
"disconnected" one 2 . 

It is important to note that in writing Eqs. 1)2. 6)1 - 1)2.8)) we have assumed that poles arise 
only from exact zero modes, i.e. that taking the limit m — * and performing the average over 
the full space of configurations commute. At fixed volume the only potential danger arises 
from the average distribution of eigenvalues near zero; our assumption holds if the density of 
eigenvalues vanishes at fixed non-zero index. Intuitively one expects that distributions of non- 
zero eigenvalues at non-trivial topology are depleted near zero. In xPT, as well as in random 
matrix theory ([H] and references therein), the densities behave as p v (A) ~ A( 2 l !y l +1+2Ar /), and 
no contribution from the non-zero modes is thus expected in the observables we consider. 

Since the zero modes are eigenfunctions of 75, the scalar and the pseudoscalar correlators 
contain the same information, 

lim (ml/) 2 (S 1 \x)S 3 (y)) = - lim (mVfC^ \x - y) , (2.9) 

to^O \ / v m— >0 

1 Since the space of lattice gauge fields is connected, different choices of D possibly lead to different assign- 
ments of index for a given configuration. 

2 The terms in Eqs. 112. 71 . 12.811 could also be interpreted as classical scattering amplitudes for pairs of zero 
modes. 
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and hence we only consider the latter in the following. Finally we note that as a consequence 
of the exact chiral symmetry maintained by the GW lattice regularization, the mass does not 
require additive renormalization and the products raP 1 need no renormalization at all. 



2.2. Chiral perturbation theory 



At large distances, the two-point correlator of the pseudoscalar density can be described by 
chiral perturbation theory. The leading order chiral Lagrangian reads 



C 



F 2 



-Tr 



Tr 



(2.10) 



where U £ SU(iVf), and 9 is the vacuum angle. This Lagrangian contains only two parameters, 
the pseudoscalar decay constant -Fand the chiral condensate E, while none of the higher order 
Li coefficients of Gasser and Leutwyler appear at the next-to-leading non-trivial order in the 
e-regime, mEV< 1. The chiral theory operator corresponding to V 1 in Eq. (|2.2j) reads, at 
leading order, 



i— Tr 
2 



T i( e io/ NfU _ u i e -ie/N^ 



(2.11) 



The correlators computed in are referred to with the notation 

Cl\x-y) ee (p I (x)P I {y)) v , (2.12) 

where I is not summed over, and the expectation value is taken at the topological index v. 

The correlators Cl I (x) have been computed by Hansen in the e-regime without fixing the 
topology [3], up to relative order 0(e 8 ), according to our conventions for the counting rules 
of the e-expansion as they are specified in Appendix A Our goal in this section is to repeat 
this calculation but at fixed topology. 

Following the notation of [Ij, the general structure of the correlator is (before volume 
averaging) , 



Ci\x) = & + ai G(x) + 0! \G(x)} + 7/ / d 4 y G(x - y)G(y) + e T 5^ ( 



where 



G(x) 



V 4±. 



1 



J n,0 



p = 2ir 



/no n 
V T ' L 



(2.13) 



(2.14) 



In dimensional regularization, G(0) = — fii/vV, with j3\ a dimensionless numerical coefficient 
depending on the geometry of the box. According to Eq. (|2.13|) . the result factorizes to terms 
representing space-time dependence, and to the coefficients Cj, aj, /?/, 7/, ej, which turn out 
to contain integrals over the zero-mode Goldstone manifold. While all these quantities depend 
on the leading-order low-energy couplings F and E, the constant Cj and the contact term e/ 
also depend on a combination of the Lj coefficients of Gasser and Leutwyler at the NNLO at 
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which we are working 4 . To avoid the dependence on these additional couplings we will only 
consider the time variation of the correlators at non-zero time separations. For convenience, 
we also average the correlators over the spatial volume L 3 . 

The various time dependences remaining after integration over the spatial volume are listed 
in Appendix lB.il the emerging zero- momentum mode integrals in Appendix IB. 2\ and the ex- 
pressions for the coefficients Cj, aj, (3i,7i in terms of the zero- mode integrals in Appendix lB.31 
For Ci the expressions are at NLO only, for the aforementioned reason. 

Taking the volume average and considering the time derivatives of the residues of the 1 / m 2 
poles, we define 



hm(my) 2 A|d 3 xCr(x) = l -A'{t), 

A c 

lim {mVf— / d 3 xC°°(x) = N v A'(t) + N 2 A'(t) . 

m^O dt J 



(2.15) 
(2.16) 



The spectral representation of Eq. (|2.6j) and the definitions A(t) = fd 3 x.A u (x), A(t) 
Jd 3 x A u ( x) , then imply that at large t, 



A'(t) = A'(t), A'(t) = A'(t) . 



(2.17) 



These constitute our basic relations between the zero-mode amplitudes and the pion decay 
constant in the chiral limit, F, once we spell out the right-hand sides. The latter actually 



vanish at the lowest order where the undifferentiated quantities A{t) = |z^|L 3 ,^4(t) 



2 L 3 



are constant, matching the sum over volume of the zero-mode expressions in Eqs. (|2.7|) . ()2.8|) 
3 . On the other hand, given the expressions in Appendix B we obtain, at NNLO, 



F2A 'V = |{(l + iVfHKW + ^(r)} 

T 2 



2\u 



FA ' (t) = -m (N ' +M ^ (r)+ N^v 

where r = t/T, and the functions appearing are given by 



Hot) 



(2.18) 
(2.19) 



H 2 (t) 



H2(r) 



+ 



(i + jv f M)jV( 
i 

2 



h'^r) + N{(6 - Af )M + 4 + Nf(2v 2 - 1) h' 2 (r) 



N{(2 - 1 N?)\v\ + 1 + -A^ f 2 



-(AT f + luDNf^^hiir) + [(4 + iV f 2 )|z,| + 2AT f (2 + v 2 ) 



h' 2 (r) 



+ 



'l + \N 2 )\u\+2N i - l -Nf]g' l (r) 



(2.20) 



(2.21) 



The functions hx,h 2 ,gi are defined in Eqs. (|B.5|) (|F3.7|) . Note that only the low-energy cou- 
pling F appears here. A non-trivial check of these formulae is that, for \v\ = 1, they satisfy 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
t/T t/T 



N f =3 N f =3 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



t/T t/T 

Figure 1: The NLO (solid, volume-independent) and NNLO (dashed for L = T = 2 fm, 
dotted for L = T = 3 fm) predictions for F 2 A'(t), F 2 A'(t) at \u\ = 1,2 (thick, thin), for 
N { = 2 (top) and N f = 3 (bottom). We have chosen F = 93 MeV. 

-A'(t) + -A'(i) = for any Nf, as must be the case since the sums in Eqs. (|2.7(l and ()2.8(l are 
identical if there is only one zero mode. 

Fig. H] shows the NLO and NNLO results for F 2 A'(t), F 2 A'(t) as a function of time, for 

3 Provided that the probability of having zero modes of both chiralities is zero. 



6 



iVf = 2 and iVf = 3 and two volumes. Considering, say, the slope of the curves at around 
t/T = 0.5, the NNLO correction is ~ 50% of the NLO term in the smaller volume shown, if 
\v\ is not too large, and then decreases in larger volumes as ~ 1/y/V. 



3. Quenched lattice determination of the low-energy couplings 

In this section, we move on from the full theory, which at present is not easily accessible to 
lattice techniques, to consider its quenched approximation. We derive the quenched chiral 
perturbation theory |231 124j predictions for the pseudoscalar correlation functions of zero- 
mode eigenfunctions and compare them with numerical results obtained in quenched QCD 
with the overlap Dirac operator. 



3.1. Correlators in quenched chiral perturbation theory 



The predictions obtained with xPT, Eqs. (|2.18j) and (|2.19|) . diverge in the formal limit Nf — > 0. 
This indicates that the results will be substantially modified in the quenched theory. Our 
working hypothesis is that correlators of the form of Eq. (|2.12|) can nevertheless, at large 
distances and in a certain kinematical range, still be described by an effective chiral theory, 
called quenched chiral perturbation theory (QxPT). 

The most important difference between QxPT and xPT is that the singlet field <l?o ~ 
IndetC/ cannot be integrated out in QxPT |231 124j . The corresponding chiral Lagrangian 
may then contain all possible couplings of the singlet field and the theory loses much of its 
predictive power, unless an additional expansion in 1/N C is carried out. In this case, the 
analysis of the relevant operators follows very closely the analysis of the generalized chiral 
theory, including the 7/ in full QCD |2£ 



|27j . and is reviewed in Appendix A The presence 
of new couplings implies that e-counting rules have to be defined for them. There are several 



possibilities, as we also discuss in | Appendix A[ We choose one that has not been considered 
previously, to our knowledge, for reasons that will presently become clear. 

In the so-called supersymmetric formulation, the quenched chiral Lagrangian at the order 
we are working reads 



QxPT 



F 2 

— Str 
4 



d^ud^u- 1 



UgU 



U- X U2 



-Str 



UgU + U~ X U7 X 



+ 



2N r 



a 
2N, 



(3.1) 



where U £ G1(JV V |JV V ) [2H1, Str denotes the supertrace, $0 = "f"Str [— iln(U)} and the vacuum 
angle 8 appears as Ug = ex.p(i6Ipf v /N v ), where In v is now the identity in the physical N v x A^ v 
"fermion-fermion" block and zero otherwise. The matrix Ug commutes with all the flavour 
group generators T 1 , which are also assumed to be extended to become 2N V x 2A^ V matrices, 
with only the physical block non-trivial. Besides F, E, the quenched Lagrangian in Eq. (|3.1|) 
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contains now three additional parameters: K, tuq/Nc and a/N c . At the same order as 
Eq. (|3.1j) . the operator corresponding to Eq. (12.21) becomes 



i-Stx 



T'lUeU 



i^oStr 



T T (UgU + U- 1 U- 1 



(3.2) 



In the e-counting we have adopted, the mass parameter related to the sing let field, ml/N c , 
will be treated as a small quantity of 0(e 4 ), so that only the first order in it needs to be 
accounted for. The reason is that this guarantees that the non-zero mode Gaussian integrals 
over the graded group, performed according to Zirnbauer's prescription |28M29| . are formally 
well defined. This counting also automatically implies that 



(4irF)< 



(4vr) 2 iV c 



nir. 



(3.3) 



which is the window where QxPT should converge. Indeed, quenched corrections increase in 
size with the volume in contrast with the unquenched case where they decrease: contributions 
of the form m^y/V fN c ~ {v 2 )/\/VF 2 become large if we do not satisfy Eq. (|3,3[) . In the real 
world, obviously, 1/N C is not tunable, and a phenomenological justification for the counting 
introduced is simply that it seems to be able to describe our data, as shown in the next 
sections. 

Correlators of the form of Eq. I|2.12[) again factorize into two types of pieces, space-time 
integrals and zero- mode integrals. In the quenched theory the zero-mode integrals can only 
have terms oc N v (from the connected contraction) and oc N% (from the disconnected one). 
The connected contraction then directly determines the result for the non-singlet correlator. 
The two parts can be determined as discussed in |30j : the former by using the replica for- 
mulation [35 and the U(iVf) integrals that already appeared in the full theory, the latter by 
carrying out the full computation of the zero-mode integrals for iV v = 1, and subtracting the 
connected part. Therefore, it is enough to consider C°°, for a general N v , and deduce C^ a 
from the part oc N v in C®°. 

Generalizing Eq. ()2.13() . the overall form of the answer now is 



C™(x) = C + a G(x) + a' E(x) + (3 G{x) + p' G{x)E(x) + ffi E(x) 



(3.4) 



+ 



d 4 y \ l0 G(x - y)G(y) + j Q G{x - y)E{y) + ^E(x - y)E(y)] + e ^ 4 )( 



where, instead of E(x) = G(x)/N{ as in the unquenched theory, we now have 



E(x) = 

Here G{x) is defined in Eq. 1)2. 14j) . and 



W G(x) + ^F(x) 



1 / , „ \ < => V P' X 



(3.5) 



(3.6) 
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The additional time-dependent functions appearing in the quenched case, owing to the 
function F(x) in Eq. ()3.5|) . are listed in Appendix lC.il The quenched zero- mode integrals 
are discussed in Appendix IC.2| and the expressions for the coefficients in Eq. 1)3.4 J) , in terms 
of the zero-mode integrals, in Appendix IC .31 

Collecting everything together, we obtain for the objects in Eqs. ()2.15)) . ()2.16j) . 



F 2 A'{t) = 2\v\ 



\u\h' l {T) + 



a 



2KF 



P 



+ 



T 2 



F 2 A'(t) 



+- 



F 2 V 
2\u\ 

rp2 



2v 2 + 



\2iV c S 
;-2(z, 2 ))/ i ' 2 (r) + 



F 2 VV 
T 2 



h' 1 {r) + \u\ 
— \v\ 



a 
2N~ r 



2KF 



2F 2 V 
Pi 



2\v\(v 2 ))ti 2 (T) + \v 



F 2 VV 

rp2 



F 2 V\ 3 

We have used here the Witten-Veneziano relation ijiqF 2 



h[(r) 

g'xir) 

h[(r) 
9i(r) 



2F 2 V 
4N c (u 2 )/V 



(3.7) 



(3. 



which is exact at 



1, 



this order, where (z/ 2 ) /V is the topological susceptibility. It may be noted that for \i> 
A' it) + A'(t) = 0, as should be the case. We observe that there are three independent low- 
energy parameters entering the expressions: (i^ 2 ), the combination a/2N c — 2KF/T,, and 
F; we thus set, without loss of generality, K = 0. Obviously a simultaneous determination 
of three parameters from the zero-mode eigenfunctions will be more difficult than in the 
unquenched case, where only F appears at this order. 

In Fig.Hwe show the NLO and NNLO predictions for F 2 A'(t) and F 2 A'(t), Eqs. (J377J) and 
()3.8|) . for L = 1.6 and 2.0 fm. Considering, say, the slope of the curves at around t/T = 0.5, 
the NNLO correction grows to ~ 50% of the NLO term in the larger volume. 



3.2. Simulation details 

We have performed a lattice simulation in the quenched approximation, using the overlap 
Dirac operator for the fermions The topological index and the zero-mode eigenfunctions 
are computed as proposed in ^7] on thermalized configurations for two physical volumes 
and various lattice spacings. Only sectors with topology \u\ = 1,2 are considered. Tabled 
summarizes the simulation parameters; the same configurations have previously been analysed 
in a different context ,8,. 

From the zero-mode eigenfunctions, we compute the volume average of the correlators in 
Eqs. 1)2.7)) and 1)2.8)1 . There is a good signal in all illustrated in Fig. EJ QxPT 

predicts that at non-zero times these correlators should behave like polynomials in time. We 
thus consider a Taylor expansion around the mid-point, r = 1/2. Denoting z = r — 1/2, we 
define the coefficients D v and C v as 

^A'(t) = D v z + C v z z + 0(z 5 ) , (3.9) 
72-4' (t) = z + C u z 3 + 0(z 5 ) . (3.10) 
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Figure 2: The NLO (solid, volume-independent) and NNLO (dashed for L = T = 1.6 fm, 
dotted for L = T = 2.0 fm) predictions from Eqs. (|3.7j) . (|3.8j) . The QxPT parameters are 
a/7V c = 0, F= 115 MeV, (f 2 )/V = (200 MeV) 4 , and \u\ = 1,2 (thick, thin). 



Lattice L/a r /a L[fm] N meas (\v\ = 1) iV meas (|^| = 2) 
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5.8458 


12 


4.026 


1.49 


880 


696 


Bi 


6.0 


16 
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1.49 
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B 2 


6.1366 
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1.49 
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213 


Co 


5.8784 


16 


4.294 


1.86 


229 


186 


Ci 


6.0 


20 


5.368 


1.86 


83 


78 



Table 1: The simulation parameters (cf. ref. [B]). Here a is the lattice spacing, r$ is the 
Sommer scale pS; . tq = 0.5 fm, and iV meas is the number of configurations. All lattices are 
symmetric, T = L. 



With a simple linear fit we can then extract the parameters D v and C u on jackknifed config- 
urations. Table[2]shows the results of these fits in the time interval (t m inj T/2). The data are 
modelled very well by the fits, and also the dependence on the choice of £ m j n is insignificant. 

It is clear from Table |2] that only the D u coefficients can be extracted from the data in a 
reliable way. The errors on the C u coefficients are large and their central values vary quite 
significantly with the lattice spacing. This is to be expected since the C v coefficients are 
more relevant at short distances and so will also be more sensitive to cutoff effects. For this 
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Figure 3: The numerical data corresponding to A'(t)/L 2 , A'(t)/L 2 , from the B lattices. 
Where not visible the statistical errors are smaller than the symbols. The left plot is for 
\v\ = 1, the right one for \v\ = 2. Comparing the slope at t/T = 0.5 with the xPT predictions 
shown in Fig. allows us to estimate (FL) 2 . 
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10.9(2) 


- 9.1(4) 


- 7(1) 


-24(1) 


19(3) 


Bi 


4a 


0.219 


3.9(2) 


10.8(4) 


-11.1(9) 


- 5(2) 


-11(3) 


26(8) 


B 2 


•5a 


0.225 


3.8(3) 


10.7(4) 


- 9.7(9) 


- 3(3) 


-13(4) 


8(8) 


Co 


4a 


0.219 


3.4(3) 


8.9(5) 


- 9.2(9) 


2(3) 


- 4(5) 


16(9) 


Ci 


5a 


0.225 


3.9(5) 


9.4(7) 


- 9.1(16) 


-11(5) 


7(7) 


2(19) 



Table 2: The D u and C u coefficients from the fit. The numbers in parentheses indicate the 
error of the last digit. 



reason, we restrict ourselves to the D v coefficients in the following. 



(|3~7|) and (pTS]l and a matching with Eqs. (|33|) 
2KF fa 



3.3. Analysis of the data 

A Taylor expansion of the functions in Eqs. 
and (|3.10|) gives 

D u = TT ± T L T {\v\ + 



(FL) 2 V 1 2N C 
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Lattice 


dof 


x 2 - 

A.min 


FL (B) 


a/N c 


X 


B 


6 


5.8 


(0.84,0.92) 


(0.3,1.2) 


(0.04,0.08) 


C 


3 


2.3 


(0.76,1.01) 


(-0.6,2.6) 


(0.04,0.08) 


B+C 


12 


8.8 


(0.83,0.90) 


(0.3,0.9) 


(0.05,0.08) 



Table 3: Results from the global fits. The intervals are the projections of the 68% confidence 
level contours. 



+ 



T 2 



3 J 2 

£>„ = / nr \Q ^ -1 - 



(FL) 2 l 1 ] v 2iV c S F 2 ^ 

T 2 



— \v\ 



F 2 F J ' 



■1 



(3.11) 



(3.12) 



F 2 F J ' 

where we have written h' 2 (T) = + 0(z 3 ), (^(t) = 712 + 0(z 3 ), and 

C 2 = 7 , 71 = -— + -V o~ w -0.0571276522 , for T = L , (3.13) 

24' 12 2^ sinh 2 (|p|/2) > K ) 

with |p|=2vrTE? =1 n 2 ] 1 / 2 /L. 

The quantity (u 2 ) in Eqs. (|3.11|) and (|3.12|) has recently been computed with high accu- 
racy [SI . In ref . [S] the results obtained at several lattice spacings were consistent with a 
well-defined continuum limit [SHI, giving \ = r o(^ 2 )/^ = 0.059(5), where tq = 0.5 fm [3*2*] . 
In order to determine the other parameters F,a/N c , we need to fit for them simultaneously, 
but also take into account the error in the determination of x- By comparing the results for 
the D u coefficients on the different lattices, cutoff effects are seen to be negligible within the 
statistical uncertainty. For this reason we do not attempt a continuum extrapolation here 
and simply consider the data at different lattice spacings as statistically independent. Since, 
on the other hand, the value of x cited above is the result of a continuum extrapolation, we 
will assign new error bars to it, large enough to also incorporate the finite lattice spacing 
values from [H]: xo = 0.059(10). We then perform a x 2 minimization in the three-parameter 
space (F,a/N c ,x), with x added to the x 2 function as 8x 2 = [{x ~ Xo)/^Xo] 2 - 

We have performed three fits: the B lattices, the C lattices and their combination, taking 
into account the ratio of physical volumes, L^/L^ = 16/20. The values of Xmin anc ^ ^ ne 
projections of the 68% confidence level contours onto the different parameter axes can be 
found in Table 01 The quality of the fits is good, with Xmi n /dof < 1.0 in all cases, and the B 
and C lattices give rather compatible results. 

It is interesting to contrast this situation with what it would be in full QCD, where the 
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1.0 1.5 2.0 2.5 3.0 

L [fm] 



Figure 4: The coefficients D±, D2 and D2 as a function of the box size for Nf = 2 (solid), 
TVf = 3 (dashed) and quenched (dotted). The quenched parameters are chosen at their best 
fit values according to Table while F = 93 MeV for Nf = 2, 3. The symbols denote the 
averages of the data points from the B (L( B ) = 1.49 fm) and C (L^ = 1.86 fm) lattices. 



D v coefficients only depend on the decay constant F: 



D„ 



+J^L;i[l + \v\Nf)(l-Nf h 



Nf(FL) 2 



+ C 2 ( (6-^)^1 + — + Nf{2v 2 -\ 
2\v\ 



+ 



T 2 



Nf(FL) 2 { 1 "V F 2 \fV ) F 2 V 



(N { + H)( 1 -Nf 
+ v>( I + 2//- - Nf + ^\v\ + Nf\u 



F 2 VVJ F 2 V 
1 

01 



71 



2 + Nf 4 - Nf 
2Nf + 2 ' 



(3.14) 



+ 



f4-Nf t 2 + Nf 
2Nf ' 



n 2 



+ 



(3.15) 



These expressions show reasonable convergence (in the sense that the NNLO correction is less 
than 50% of the NLO term) only at L > 2.0 fm for F = 93 MeV, and in order to push the size 
of the correction below 30%, one would need to go to L > 2.5 fm. In this case we might expect 
a systematic uncertainty in the determination of F of about 5%, and statistical uncertainties 
would reach the same level if D u could be determined using ~ 100 configurations with \u\ = 1. 

The full theory formulae, Eqs. (|3.14j) . (|3.15|) . possess the feature that the NNLO corrections 
come with negative relative signs, such that the expressions are almost independent of FL at 
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around FL ~ 1, and their absolute values have an upper bound at this order. For illustration, 
we show the full predictions for F = 93 MeV in Fig. |3J as a function of the box size (for a 
symmetric geometry, L = T). Also shown are the quenched data as well as the quenched 
predictions. Because of the mentioned near-cancellation, the full predictions at this order 
could not be moved significantly closer to the quenched data points by tuning F. In any 
case, as already mentioned, they show reasonable convergence and can thus be considered 
self-consistent predictions only for L>2.0 fm. 

4. Conclusions 

Approaching the chiral limit has remained a long-standing challenge for lattice QCD for 
many reasons, among them that finite-volume effects become large for very light pseudo- 
Goldstone bosons, and that the Dirac operator develops very small eigenvalues. It has been 
the purpose of this paper to elaborate on the fact that at least these particular problems can 
be overcome: for instance, the Dirac operator eigenfunctions associated with the exact zero 
modes encountered in gauge field configurations with a non-trivial topology at finite volume, 
can be used to extract physical information concerning the chiral limit of the infinite-volume 
theory. 

More precisely, we have shown that certain classical scattering amplitudes of the zero- 
mode eigenfunctions measured at finite volumes, Eqs. (|2.7|) . (|2.8|) . allow the extraction of the 
infinite- volume pion decay constant, via the relations in Eq. (|2.17|) . We have worked out these 
relations to NNLO, Eqs. (|2~T8|) . (|2~T9|) . finding that the convergence of chiral perturbation 
theory seems reasonable for these observables, provided the volume is above ~ (2.0 fm) 4 . 

Finally, to estimate the practical feasibility of using such relations, we have carried out 
lattice Monte Carlo simulations in the quenched approximation, using overlap fermions. We 
find a good signal for the observables in Eqs. (|2,7|) . H2.8|) . shown in Fig. |21 Matching with 
chiral perturbation theory predictions relevant to the quenched approximation (which show 
reasonable apparent convergence only in volumes between ~ (1.0 fm) 4 and ~ (2.0 fm) 4 , 
in marked contrast with the unquenched case), we find that the pion decay constant, to 
the extent that it is a well-defined quantity in this case, can be extracted with about 5% 
statistical accuracy, utilizing a few hundred configurations with non-trivial topology. The 
number we obtain in volumes ~ (1.5 fm) 4 is in the ballpark of 115 MeV. 

Our result for the pion decay constant in the chiral limit is larger than what one would 
expect in Nature: conventional xPT in infinite volume [2S3 yields F ~ 87 MeV, if the 
standard phenomenological values for the 0(p 4 ) Li coefficients of Gasser and Leutwyler are 
inserted [31]. The fact that our quenched calculation seems to overestimate F is consistent 
with other recent quenched results for the physical F n in the continuum limit, however |35|- 
[3*7| . For instance, the results of [SHI EH] imply that the quenched F n is 10% larger than 
the experimental value, if the scale is set by r$ _32.. On the other hand, these standard 
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approaches (unlike ours) have to rely on quenched chiral extrapolations in the light quark 
masses, which introduce significant systematic uncertainties of their own |38| . 

On the side of our approach, it is conceivable that a smaller value for F could be obtained 
by going to larger volumes. As we have discussed, however, the peculiarities of the quenched 
approximation imply that the volume cannot be increased too much, since the convergence 
of quenched chiral perturbation theory soon deteriorates. Therefore, a systematically im- 
provable determination of F by using our method (or any other) lies beyond the quenched 
approximation . 
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Appendix A. Large-iV c counting in the e-regime 



Large- N c counting in the context of chiral perturbation theory has been analysed in detail in 
ref. [22j. The same general discussion goes through in the full and in the quenched theories, 
with the replacements in the latter that — > U~ l ,e ie l Nt — > U$,Ti — > Str. For simplicity, 
we will mostly use the notation of the unquenched theory here, indicating then the important 
point at which differences arise between the two cases. 

In general, the chiral theory including the singlet is, at leading order in the momentum 
expansion and to all orders in 1/N C , of the form 

- {V2($o)Tr[e"' N *UM] + VJ($ )Tr [M^U^e^ Nt ] } + ^ 3 (*o) (d^ ) 2 + ... , (A.1) 

where <&q = — *f"Tr InU and M = diag(m,m, ...). 

The Lagrangian in Eq. (jA. 1|) contains an infinite number of parameters, since the potentials 
Vi($>o) are arbitrary functions, with the only constraint from parity that Vi{— <3?o) = Vi(&o)i 
for i = 0, 1, 3 and V^(— $o) = ^2*(^o)- It can be shown, however, that they involve a specific 
power series in 1/N C ([2S1-|2Z1) and references therein). Noting that the field 4>o of [2Zj * s 
4>q = 2^q/{F^/W c ) in our notation below, the structures arising are 



^o($o) 



1 m§ 



V.($n) = 




2 

a 
2N f 



K 

+ i^=$o + 



+ 





(A.2) 
(A.3) 
(A.4) 
(A.5) 



where all parameters introduced (wq, F, S, K, a) are assumed not to scale with N c . Inserting 
the specific terms shown here into Eq. (|A.1|) . one obtains the theory up to 0(1/N C ). In the 
following, we denote 



F 



F 2 N C , 



£ = T,N C , K = 



K 



(A.6) 



In order to define a formally consistent framework, it is convenient to now combine the 
momentum and 1/N C expansions. Following |27j . we may choose 

P 2 ~ ^ ~ <? ■ (A.7) 

The defining property of the e-regime is that the pions are off-shell since the momenta are 
fixed by the size of the box, p ~ 1/L ~ 0(e), while the quark mass is small, such that 



mT,V < 1 . 



(Ai 
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Given Eqs. (TOl) . (|A~7|) . we are thus led to the rule 

,6 



m ~ e 



As usual, the field configurations are factorized into zero- momentum modes Uq, $o an d non- 
zero modes £, «l>o> 



C/(x) = C^(x)C7o , ^(s) = exp 



F 



, $ (a:) = *o(a:) + $o 



(A.10) 



where Jd 4 x^(x) = and $o = Tr£. The counting rules for the non-zero modes, which are 
treated perturbatively, are 

^~p~e, $o~P~ e - (A-ll) 

For the zero mode Uq we have Uq ~ 1, while the counting of $o is to be determined presently. 

Indeed, let us consider the terms involving explicitly the flavour singlet zero mode $o- We 
are interested in carrying out the computation up to and including 0(e 4 ), and the terms 
potentially of this order, after integration over space-time, are 



/ 



d 4 x£ x PT 3 + 



1 m 



2N C 
i mKV <l Tr 



oK)o(^)on 



(A.12) 

O(^)0(l>o)0(e 3 ) . (A.13) 



Moreover, we want to carry out the computation at a fixed topology; performing the integral 
over 9 with the weight exp(i9v) introduces (after a shift) effectively one more term, 



2iu - 

d 4 x£ xPT 3 — $o ~ O($ )O(e). 
b 



(A.14) 



Once the integral over 9 is converted to a Gaussian over 

$o, Eqs. (lA~12l) . (lA~T4l) tell that the 
saddle point is at leading order in e at $o ~ N c u / {m^FV) ~ 0(e 3 )/0(mQ). Thus, we have 
fixed also the counting of <&o. 

We can now collect together the full theory at fixed topology. The factorized part of the 
zero-mode partition function becomes 



Z v (ji) oc e ^ 



^tv./jt , tt\\ I o,, ^Tv/rL rn 



t/o e U(Ar f ) 



det^C/oexp ^Tr(C/ + f/ T ) + 2 ^ >. 



where 



(v 2 ) _ mlF 2 



Tr(£/ -£/o T ) ,(A.15) 



(A.16) 



V 4iV c 

The first term in the exponent is 0(1), while the latter is, as follows from Eq. ()A.13|) with the 
given estimate of <i>0; ~ 0(-ftT)0(e 6 )/0(mQ). The non-zero momentum modes, on the other 
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hand, are described by 




(A.17) 



(A.18) 



(A.19) 



(A.20) 



To finalize the setup, one has to decide what kind of counting rules are chosen for the 
parameters rriQ,a,K. For simplicity, we will assume that 0(a) ~ 0(K) ~ 1. The counting 
of Trig then leads to three distinct possibilities: 

1. In the unquenched theory, m 2 , can be taken as "large", say m 2 , ~ e~ 2 . Parametrically, 
then, ijiq/Nc 3> p 2 . In this case the term in Eq. (|A.20|) dominates the action: the 
non-zero modes $o (representing the rf) are heavy and can be integrated out, resulting 
in the simple usual chiral theory following from Eq. ()2.10|) , 

This choice is not available in the quenched limit, however: the field $o cannot be 
integrated out [231 [2U, but has to be treated as a light degree of freedom. Therefore, 
it is convenient to assign a different counting to it. 

2. In the "standard" version of quenched chiral perturbation theory, one chooses TOq ~ 
0(1), such that mg/iV c ~ p 2 29 . Then Eq. ()A.20j) is of the same order as the usual 
kinetic terms following from Eq. ()A.17|) . The term with K in Eq. 1)A.15|) . on the other 
hand, can be neglected, since it is 0(e 6 ) with this counting. 

This standard counting suffers from some problems, however. First of all, the integrals 
over the graded group of the supersymmetric formulation do not appear to be, strictly 
speaking, well defined ^29.^, because the masses related to quadratic fluctuations, treated 
according to Zirnbauer's prescription [2Hj, are not positive-definite (cf. Eq. (3.8) in |29|). 
Second, QxPT leads to a perturbative expansion parameter ~ m^/ {N c p 2 ). With the 
standard counting this is of order unity, formally spoiling the convergence. 

3. Because of the problems of the standard counting, we will consider an "alternative" 
counting here. In the alternative counting, uiq is treated as a small quantity, say 
m 2 , ~ e 2 . Then Eq. (|A.20|) is formally a perturbation and the leading order quadratic 
form is well defined. This counting makes also explicit the fact that Q^PT should only 
work in the window of Eq. ()3.3|) . With the choice to 2 , ~ e 2 , contributions from the 
coefficient K should be kept in the results. 

In the real world, obviously, 1/N C is not tunable, and not necessarily small. Therefore, the 
success or failure of the frameworks described remains ultimately to be judged empirically, 
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by comparing them with data. The expressions below are, for generality, for the "standard 
counting", while in the actual text we only showed truncated versions, where terms of higher 
order according to the "alternative counting" had been dropped. 



Appendix B. Detailed results for full chiral perturbation theory 
B.l. Space-time integrals appearing 



After integration over the spatial volume, the time dependence of Eq. (|2.13j) appears in the 
forms 



d 3 x 



Jd 3 xG{x) 

d 3 x \g{x)] 2 
|d 3 x Jd A yG(x-y)G(y) 



Tht{T) , 
-T 3 h 2 (r) , 



where 



hi(r) 
9i{t) 



1 
2 
1 

24 



r 

2 



r 2 (r - 1) 



1 

12 
30 



Mt)] 2 + £ 

n^O 



cosh (|p|(r- 1/2)) 



i 2 



2|p|sinh(|p|/2) 



Here 



T 

IPI = 2vr z 



E 

.i=i 



»7 



1/2 



(B.l) 
(B.2) 

(B.3) 
(B.4) 

(B.5) 
(B.6) 
(B.7) 



(B.8) 



B.2. Zero-mode integrals appearing 

The zero-momentum mode integrals at fixed topology are related to the partition function 



Z„(/i) = / det^ e^ ReTtl/o , 

Ju eV(Nf) 

where [i = raTiV . The value of Z v {\i) is known |39[ I18j to be 

Z y (ii) = det[I v+ j-i((J,)], 



(B.9) 



(B.10) 
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where the determinant is taken over an Nf x Nf matrix, whose matrix element is the 
modified Bessel function I u+ j-i. We will express our results in terms of the derivatives of 
this partition function, in particular 



_ E„(/a) _ 1 d 



At small \x and non-zero v, 



<T v (pt) PS . 



Expectation values are denoted by 



<•■•>£ 



c/oeU(ATi 



f) (...)det^ e 



iu Re Tr U 



u ev(N { ) 



det^C/ e' iReTrC/o 



(B.ll) 



(B.12) 



(B.13) 



both the superscript and subscript in (...}(? are often left out. 

Given these definitions, all the emerging expectation values can be computed analytically, 
using the techniques discussed in Appendix B of The small-// (small-m) limits are then 
obtained by using Eq. (|B.12|) . We give here a complete collection of the integrals appearing, 
up to third order in the matrices Uq,Uq. The expectation values for complex-conjugated 
operators are obtained from those shown simply by v — ► —v. For the small-// limits we only 
show the values of order for n powers of Uq, U^: 



N 



v 

Ov 

— (1/ -u) , 



<Tr([/ 2 )> = Nf 



1 



2(N { + v) 



v 

<?v 



.™*.(N t + v){\u\-u) 



(TrC/oTrC/J) 



((Trt/o) 2 ) 



2, 



Ni\o f v + N t <% + — -N£ 1 

N f 
x -=■ , 



2Vf 



<t[, + A^ 2 - (1 + 2iV f i/) — + (2 + ATfi/) 



—±(l+Nfv)(\v\-v) 
A* 



(B.14) 

(B.15) 

(B.16) 

(B.17) 

(B.18) 

(B.19) 
(B.20) 

(B.21) 



Af 

+ 



^(2^ 2 + 5iV f i/ + 2 A" 2 + 2) - -(2A/f + 3i/) 



/i 



1 + ^{2v 2 + 6N f v + 2Af + i; 



<7v 



K + tyo*) (B.22) 
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(TtU Tr(U$)) 



((TrUoY 



((TrU^)(TrU f 



4Nf 

t? 

Nf 



{l + N? + ?,Niv + v 2 ){\v\-v) , 



(B.23) 



2v[ 

T? 



(4 + iV f 2 )i/ + 7V f (4 + z/ 2 ) 



2 + N { u 



Ni + Jfi (( 3 + 2N t> + N t ( 3 + 21/2 
Nf + v 



-{a' v + NfCxl) 



AN f r 



2JV f + (2 + Nf)v + iV f ^ 2 - v 



(Tr^Tr^ 2 )) = iV f 



2NfV 7 



(iVf + 1/) + 



+ 



iV f --^(JV f + 

^ 2 



2 + iVfi/ 
A* 

2(Nf + t/) 



(a^ + iVfa 2 ) 



iV f 



Nf 



:(8 + <SNfV + iV f V) + -^(1 + 3NfV + N(v*)a v 



-(1 + JV f i/)(oi + iV f a 2 ) + ^ + 3^^ + iV f V 3 



(2 + 3iV f i/ + iV f V)(M -i/) 



4Nf 

r I/ 2 1 
JV f iV f — (2 + Nfu) -—(1 + N t v + Nfv 2 )a v 
L fx 6 H A 

+ 1 ~ NiU K + ^) + < + 3JV f ot^ + Nfa'l 
Ox — . 



(B.24) 
(B.25) 

(B.26) 
(B.27) 

(B.28) 
(B.29) 

(B.30) 
(B.31) 



It may be noted from the small-// expressions that plenty of degeneracies emerge if we put 
Nf — > 1: this is simply because taking a trace has then no meaning. 



B.3. Results for the coefficients in Eq. (|2.13|) 

Let us define 



jV f 2 -lG(0) 
iV f F 2 



1 + 



Af - 1 ft 



iv f F 2 ^/y. 



and 



For the coefficient Cj defined in Eq. (|2.13j) . we then obtain (at NLO) 

2\ fjf 



(B.32) 
(B.33) 

(B.34) 
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where / is not summed over. For 1 = (flavour singlet), the result is immediately related to 
the expectation values listed in Eqs. (|B.14|) - ()B.31|) ; for I = a, one can make the connection 
by using independence of a and the completeness relation 

1 1 1 

( T i a j T kl F ijkl) = N 2 _ 1 J2 ( T tj T kl F ijkl) = 2(jV| - 1) ~ W { Fiijj ) ' ^ B ' 35 ^ 



We show explicitly only the small-// limits here: 

T, 2 Nf\u\ 



Cn 



C a 



A* 

T?\v\ 

"V" 



(1 - N t \v\ 



For aj, we obtain, in a similar way, 



Oil 



2F 2 
x< Tr 



_ N 2 - 2 G(Q) 
iVf F 2 " 

(U T T + T 7 ^) 2 ] - -^{Tr [T\U + C/J)] }' 



2\ A* 



For small /i, 



2S>| 
F 2 /x 2 

F 2 fi 2 N { 



1 + Ni 



1 + N { 



Nf 

F 2 
G(0) 



F 2 



(1 - iv f 2 ) , 
(i + JVfH) 



For Pi, we obtain 



N? 



4/V f 



-Tr 



lTT l? 



For small //, 



Po 

Pa 



T?\u\ 

FV^f 
Y?\v\ 
2F 4 fi 2 N 2 



(N 2 - 1)(N 2 



(UqT 1 - T I U { 



2NM) , 



iY f 2 + 2-7Vf(iV f 2 -4)H 



Finally, 77 reads 
li 



Y? I f r , - 
2F A V . 

- M Tr f(C/ + ulWoT 1 + T'C/J) 2 



:/ {Tr [T J (U + l/J)] } 2 - iV f Tr [(U T J + T 7 ^) 



2/i 
'iVf 
Aj 

iV 2 



Tr 



(C/ + c4)(^oT J + r J ?7o)l Tr h 1 (U + tf ) 



Tt(U + uI){Ti [r J (^o + ^ + )]} : 



2\ V 



(B.36) 
(B.37) 



(B.38) 

(B.39) 
(B.40) 

(B.41) 

(B.42) 
(B.43) 



(B.44) 
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For small n, 



70 
7a 



2E 2 H 
F 4 Vn 2 N { 
X 2 M 



(iV f 2 - l)(iV 2 



iV f 2 (l - 2z/ 2 ) 



-4-2JV f M) , 
4 + iV f (iV 2 -6)\v\ 



(B.45) 
(B.46) 



Appendix C. Detailed results for quenched chiral perturbation theory 
C.l. Additional space-time integrals in the quenched theory 

Apart from the integrals in Eqs. (|B. 1|) (jB .4|) . Eqs. Q3.4|) and (j3.5|) imply that in the quenched 
case we need, in general, the following further ones: 



y*d 3 xF(x) 
Jd 3 xG(x)F(x) 



d 3 x 



F(x) 



-T 3 h 2 (r) 
T 5 h 3 (r) , 



Jd 3 xJd 4 yG(x-y)F(y) 
J d 3 x J d 4 y F(x - y)F(y) = -T 7 h A (r) 

where we have defined (following the notation in Eqs. ()B.5|) - (|B,8|) 1 
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C.2. Quenched zero-mode integrals for the flavour singlets 

As discussed in the main text, in the quenched case, the results for the flavour non-singlet 
follow from those for the flavour singlet. Therefore, we only need to address the zero- mode 
integrals arising for the flavour singlets, and do not present a similarly exhaustive list as 
in Appendix IB. 21 



23 



The flavour singlets contain two parts, a connected contraction (oc N v ) and a disconnected 
one (oc Ny). The replica trick (cf. 30 ) allows to obtain the result for the connected contrac- 
tion from a certain limit of U(iVf ) integrals, discussed in Appendix IB .21 For the disconnected 
contraction, on the other hand, the zero-mode integrals have to be honestly carried out, for 
Ny = 1, using the supersymmetric formulation of QxPT. (The non-zero momentum modes 
of the Goldstone bosons can still be treated with the replica formulation |31| . and only the 
remaining zero-momentum mode integrals need to be transformed to the supersymmetric 
ones.) We first list the supersymmetric integrals for N v = 1, and then the generalizations to 
any N v obtained with the replica trick. 

Let us start with some notation. We introduce a projection operator P v , 



(Pv)ij — 



5ij, for i,j = physical flavours in the valence block, 
0, otherwise. 



(CIO) 



Using again the scaling variable fj, = mSF, all mass dependence of the results can be expressed 
in terms the same zero-mode integral as appears in the quark condensate obtained with 
Gl(l|l) 00]: 



2K 



lMK v {p) + I U+1 (/i)i^-i(M) +-, (Cll) 



where I u , K u are modified Bessel functions. We recall that, for v ^ 0, 



as in Eq. (|B.12|) . Note that, in contrast to Eq. (|C.11|) . 

(Str^o + C/Q- 1 )) = 0, 



(C.12) 



(C.13) 



and also that (Str (U - Uq 1 )) = 0. 

The zero-mode integrals for N v = 1 can be derived following the techniques discussed 
in [22], particularly the explicit parametrization of Gl(l|l). The integrals needed, and their 
small-/i limits, read (U = Uq here) 
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Using these integrals together with the iVf — > limits of the corresponding U(iVf) integrals 
from Appendix IB.2I (obtained, in each case, with the replacements Str — ► Tr , Uq 1 — ► C/q, 
P v — ► 1), we can deduce that 
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In Eq. (|(J.44|) we used Eq. ()C.13|) together with the fact that, in general, 



Str^o + f/o" 1 )]^) = (Str(C/ + ^o _1 ))( M ) + 2 ^( M } ■ 
C.3. Results for the coefficients in Eq. (|3.4|) 

Given these building blocks, we can collect our results together. In analogy with Eqs. 
I)B.33|) . we define 
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and 



mYlV . 



Then the results for the coefficients in Eq. (|3.4|) . together with the small-/i limits, read: 
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